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CTN , We propose a kinetic model of transient nonequilibrium phenomena in metals exposed to ultra- 

Q j ' short laser pulses when heated electrons affect the lattice through direct electron-phonon interaction. 

This model describes the destruction of a metal under intense laser pumping. We derive the sys- 
tem of equations for the metal, which consists of hot electrons and a cold lattice. Hot electrons 
,^^ ' are described with the help of the Boltzmann equation and equation of thermoconductivity. We 

^5^ , use the equations of motion for lattice displacements with the electron force included. The lattice 

deformation is estimated immediately after the laser pulse up to the time of electron temperature 
relaxation. An estimate shows that the ablation regime can be achieved. 
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The first theoretical prediction of transient laser- induced nonequilibrium electron temperature phenomena in metals 

I ■ was made more than twenty years agoEl. It was shown that an ultrashort laser pulse (~ 10~^^ — 10~^^ s) produces 

C^ , a nonequilibrium state of the electron gas near a metal surface. However, experimental picosecond (~ 10~^^ s) 

"t^ ' laser studies of thermally assistecLmultiphonon photoemission were unable to measure, and even failed to observe, 

I t [ this nonequilibriuni, electron stateB. This failure had a simple explanation in terms of the theory of electron-lattice 

Cd . thermal relaxation,El which yields a relaxation time of the order of Te-i ~ 10^^^ s. It was necessary to use power pulses 

S ■ shorter than t^-i. Such measurements with-iSubpicosecond (~ 10~^'p|S)r^ulses revealed a transienWiopequilibrium 

I . regime in transmittivity and IR reflectiorJa -O, giant electron emissionB -llil and the emission of lightll3 -EJ. 

^ ' We briefly summarize the physical process. The ultrashort laser pulse (At ^ lO"^** — lO^^'^ s) absorbed in a metal 

^ , raises the electron temperature T^ considerably higher than the lattice temperature because of the difference in their 

fj ■ specific heats (ce <C ce). Subsequent electron cooling results mainly from two processes, namely electron-lattice ther- 

I— I [ mal relaxation and electron thermoconductivity. These are usually modeled with a set of coupled thermoconductivity 

equations for the electron and lattice components. These equations are nonlinear and can generally be solved nu- 

K^ ' merically, yielding the electron temperature relaxation. The solution also shows that the subsequent ablation regime 

j«^ , can be achieved, which involves the "cold" destruction of a metal into the parts consisting of different phases. "Hot" 

^^ ■ destruction, namely melting, can also be studied with the help of this solutiontj. However, such an approach has 

f^ . several shorcomings. First, the question remains as to whether the equations of thermoconductivity are still hold at 

OQ such high frequencies {^ 1/Ai). Second and more importantly, these equations can only describe the temperature 

^P ^ dynamics of a metal but not electron transport, lattice deformation, thermionic emission, etc. It is evident that a 

strict kinetic approach is needed, to describe the various transport phenomena properly and derive thoroughly the 

equation of thermoconductivityllj. 

ti . In this paper we present a theory of transient nonequilibrium phenomena in metals subject to ultrashort laser 

pulses. Our theory is based on the Boltzmann equation for the nonequilibrium electronic partition function. We 

focus mainly on times shorter than the electron-lattice relaxation time Te-i- Electrons therefore affect the lattice via 

'"^ ' direct electron-phonon interactions. To consider lattice deformations, we use the equations of the so-called dynamical 

^ ' theory of elasticity. Lattice deformation is due to the nonequilibrium electron state and results from the effective 

O ' "driving" force (proportional to VT^) on the lattice. This force also governs the renormalization (depending on 

Te) of the lattice constants (sound velocity and optical phonon gap). We show that the driving force leads to large 

lattice deformations, and can distroy the crystal. Tiicse results are in agreement with measurements of time-resolved 

X-ray diffraction synchronized with laser pumpingtil. A nonstationary increase in lattice parameters of Au(lll) and 

Pt(lll) single crystals was detected. Measurements of the shift and intensity variation of Bragg peaks, as well as the 

C^ ' Debye-Waller factor, anables one to separate the effects of lattice deformation and heating. The transformation of 

elastic into plastic deformation was also observed. 

The plan of the paper is as follows. In Sec. 2 we present the kinetic theory of the process under study: the 
Boltzmann equation for an electron gas and the elastic equation for the lattice are derived, along with the equation 
for thermoconductivity. In Sec. 3, the solutions of the proposed equations are found for the times of interest. 
The lattice deformation is calculated. In Sec. 4 the solutions are analyzed. The lattice deformation is estimated 
analytically in various limiting cases. The possibility of crystal destruction under laser pumping in discussed. 



II. THEORETICAL FRAMEWORK 

Let us briefly recapitulate the main equatinns of our problem. For the lattice deformation we use the so-called 
equation of dynamical theory of elasticitylla,t^ 

where p is the lattice density, Xijim is the tensor of elastic constants, and the driving force describes the effect of free 
carriers on the lattice 



^ d f 2(fip , , , . , - 
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The deformation potential Ay (p) yields the change in the local electron spectrum, 

fc(p, r, t) = \j{^)u^j{v,t). 

To find the electron distribution function /p(r, i), we use the Boltzmann equation with the electron-phonon collision 
integral 

st/p - E / (03 -S^'5(^p+fc - ^p + -i"') ((1 - /p)/p'^l"^ - /p(i - /p')(i + ^i"^)) 

with the probability of a scattering process involving a phonon of the n-th branch 

efU,,(p)fc 






pk (n) 

where e,-" and w^." are the polarization and spectrum of phonons of the n-th branch, respectively. 

Since the phonon-phonon relaxation time is large (^ 10^^^ s) compared with the times of interest, the phonon 
distribution function NJ^"' takes its equilibrium value at the lattice temperature T^, 



N^^\Tf) 



exp(t^(")/T£)- 1' 



The electron-electron relaxation time due to scattering on phonons r ~ T^^ ^ 10^'^'' s (see below) is much less than 
the characteristic time of laser pumping. Therefore the electron gas is nearly in thermal equilibrium at the local 
temperature Tg (r, t) . We seek a solution of the Boltzmann equation in the form 

fp = f^{-f^)+Xp^, (4) 

where /o is the local equilibrium Fermi-Dirac partition function and Xp is the noncquilibrium part. We obtain for 
the collision integral (||) 

.-1 (,. (Ap)\ 5/o 



St/p = St/o-r-Mxp-^j^, (5) 

where the scattering rate 

n 

The latter estimate is valid when the ion temperature T^ is considerably higher than the Debye temperature; the 
dimensionless electron-phonon coupling constant g ^ X/sp ^ 1. The brackets denote integration over the Fermi 
surface 

2dSF 

V3 '-■■■•'■ 



w(27r)'^ 



The first term in dq) comes from tfie contribution of the local equilibrium partition function: 

X [Sie, ~ e,, - 4")) + 5{e, - e,, + J^^)) . (6) 

This term describes the energy flow from electrons to phonons when they are at different temperatures. This term is 
absent if the temperatures of the electron and lattice subsystems coincide. 

The nonequilibrium part of the electron distribution function has to satisfy two conditions. The first is indeed the 
conservation law of the number of carriers: 

d^'p dfo _ 
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This expression determines the chemical potential and results in the renormalization of the deformation potential: 
A(p)-A(p)-(A(p))/(l). 
The second condition 

(2^(£.-m)Xp^=0 (7) 

enables us to define the local temperature Te (see Ref.E3), i.e., to write the equation of thermoconductivity. 
Substituting Eq. (q) into the Boltzmann equation, we get 

<^,,|,,x,Ha>M._„,_.,,p,^,£^(^,.^),sv„/f. («) 

To obtain the equation for the local temperature Te{r, t), we multiply the Boltzmann equation (@) by {Sp — fi)dfo/de 
and integrate over p. With the help of Eq. (0) we find the equation of thermoconductivity, 

dT 

Ce{Te)-^+dwq=Q-a{Te~Te), (9) 

where Ce{Te) = ■n'^{l)Tf,/i = /3Te is the electron heat capacity and q is the heat fiow: 

Using Eq. (|^), we can find the last (relaxation) term on the right-hand side of Eq. (|^). For high temperatures 
Ti,Tf. ^ ud, the electron-lattice relaxation constant a is 
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The density Q of laser energy absorbed by electrons can be taken in the form 

Q{z,t) = I{t){l-K)Ke-''\ (11) 

where R is the reflection coefficient and k is the inverse penetration depth. The function /(<) describes the pulse 
shape. 

Eqiations (|l|) and (S) must be supplemented by the proper boundary conditions. We assume the simplest geometry: 
the metal occupyies the half-space z < 0. Hence, the boundary conditions for the above equations are 



dz 



0, ^ 
dz 



0, (12) 



signifying that the heat flow through the surface and the normal stress tensor component vanish on the surface. 

We also need the boundary conditions for the kinetic equation (pi) at the metal surface. These boundary conditions 
depend on the type of electron reflection from the surface. We assume the specular reflection for simplicity. 



III. DYNAMICS OF ELECTRON TEMPERATURE AND LATTICE DEFORMATION 

The above equations are nonlinear and very complicated. However, it is possible to solve them in an important 
limiting case. Below we are interested in times shorter than the electron-lattice relaxation time Te^i '^ Ce{Te)/a. In 
this case the lattice temperature can be set to the initial temperature To, and the last terms in Eqs. ([g) and (|9|) can 
be omitted. 

To solve the system (|^), ^ for the half-space with the boundary condition (jl^), we use the even continuation of 
the temperature Te(r, t) and the partition function Xp, E^nd the odd continuation of Uz(r, t), into the half-space z < 0: 

re(z<0) =re(-z,0), u^(z<0) = -w^(-z,0). (13) 

For the parallel components Ux and Uy one must use the even continuation, but owing to the fact that the external heat 
( |l l| ) depends only on z, these components are vanish. In Eq. (|8|) we discard (xp)j which represents the 'in-term' in 
the collision integral. This term accounts for carrier conservation, i.e., for the isotropic channel of collision processes. 
Therefore it does not affect the heat flow and lattice driving force. 
The solution of Eq. (||) has the form 

t 
Xp(r, t)= J dt'Xp (r - v(< - t'), t') exp (-^) . (14) 

— OO 

where Xp is the right-hand side of Eq. (g), 

^P = -A.(p)^ + ^-(^^+v— j. (15) 

Substituting the solution dlj), ( [iq ) into the heat flow (nw and integrating over the energy variable according to 
(Pp = d{ep — p)dSF /v, we obtain 

q(r,t)^-^Mdi'exp(-^)v(A+v|-)T,^(r-v(i-0,0\. (16) 

The expression (nq) is linear in T^. It is convenient to introduce the new function 8(r,f) = T^{r,t) and take the 
Fourier transform with respect to space and time variables. Then Eq. ( [l6| ) yields the Fourier component of the heat 
flow: 

q(k, c.) = -^ ( ^ ^^/--i > 0(k, ^)- (17) 

D \ o; — vk -V IT ^ I 

Substituting this result into the equation of thcrmoconductivity (O) , we obtain its Fourier component 

(cj — vk)vk 

3 \ ' w — vk + IT' 



«— ( '^ + . ..,. , ,^_i ) e(k, uj) = 2«;(1 - R)I{Lo)Uik), (18) 



where I{uj) is the Fourier transform of the pulse shape I{t). The factor Q(k) describes the spatial distribution of the 
laser field {^ , 

c/(fc/ 



K?- -I- &f 



which depends only on fc^. Equation (^^ yields the temperature dynamics of metals under laser heating with the 
time and space dispersion. 

We now turn to the equation for lattice displacements (|l|). The driving force Gi{r,t) can be evaluated as in the 
derivation of Eqs. (ttq) and (|l^). Both the local equilibrium partition function and nonequilibrium part (0) contribute 
to the integral (y). Expanding the integrals over the energy variable in powers of Te/ep up to the second order, we 
obtain 

6 dsp \ w — vk + ir^^ 



G.(k,o;) = ---- '^ "^ i^ 9(k,c^). (19) 



In addition to the electron force (jT^) we also obtain the temperature-dependent renormalization of the elastic constants 
Kjirn (sound velocities) due to the interaction with electrons (electron loop in the phonon self-energy function) . The 
dominant contribution in the range of interest comes from the local equilibrium partition function: 

Kjlm -^ Kjlm — (Aij(p)A;„j(p)) ^ Jj-T i^ij iP)^lm{p)) ■ 

F 

The electron contribution to the sound velocity is second order in the electron temperature, As/ s ^ T^/e'^. 

Taking the Fourier transform of the left-hand side of Eq. (|]), one needs to keep in mind the singularity at z = 
after continuation (|l3| ) of the function u^- This singularity contributes the dS{z)/dZ'tej:ni in the second derivative 
(fuz/dz"^; such a term accounts for surface effects. The Fourier transform with respect to the coordinate z over the 
entire space gives 

- p{uj^ - s'^k^)uz{k,uj) = Gz{k,uj) + kC{uj), (20) 

where s = Xzzzz/p is the longitudinal sound velocity in the z-direction. In the last term, C(w) must be determined 
from the boundary condition dl3), and takes the form 

^, ^ ^ f dk Gz(k,uj) ,„ , 

^H--2...y--^^. (21) 

We next proceed to the electron temperature and lattice deformations represented by Eqs. (|l^) and ( |2(]| ) in various 
limiting cases. 

IV. SPATIAL VARIATION OF ELECTRON TEMPERATURE AND LATTICE DEFORMATION 



Equation (y_8|) describes the electron temperature, evolution under ultrashort laser heating of metals. This equation 
generalizes the usual thermoconductivity equationld. We are interested in the wave vector k, which is the greater of 
the inverse skin depth k (^ 10^ cm^^) and the electron diffusion length Vy/rtQ during the laser pulse to- In the usual 
experimental situation r^^ ^ 10^** s~^, kv ^ 10^'^ s~^, and the hydrodynamic regime nv <C t~^ obtains. Thus, one 
can omit the term kv in the denominator of the left-hand side of Eq. dlq) . The dominant contribution comes from 
the diffusion pole w ~ ru^fc^ <C t^^. Therefore, we can also omit oj everywhere in comparison with t~^ or kv. The 
solution of the thermoconductivity equation reads 

oo 

«<^>') - «» + ^ / (2';^ti",^-) ^"'''""""''""'"''-"'' <^^' 

— oo 

where the diffusion coefficient D = r(w^)/(l) is introduced. The constant 8o = T^ comes from the solution of 
the corresponding homogeneous equation, and represents the initial temperature. Evaluating the integral ( |22| ) with 
respect to w and k, we obtain 

t oo 

e(z, i) = eo + / dt' f dz'-^MM^e^p f-Sl^.^) . (23) 



-oo 



P^7r{t ~ t')D V ^{t-t')D 
We see immediately that the function (^3]) satisfies the boundary condition (0^). For the temperature at the surface 



z = 0, Eq. (|23|) gives 

t 



r|(0, t) = 7(2 + A / dt' g(0, t - t') e'''^*'erfc {^/^^Iw 
ttP J \ 





The electron temperature (EST) just after the pulse peaks at the surface: 



TL. - ''°% "^^ min (., (Dto)-'/') . (24) 



This result has a simple explanation. For short pulses Ky/Dto <C 1 the time dependence of the temperature corresponds 
to the local laser intensity at the observation point. In the opposite case, ny/Dta S> 1, the temperature distribution 
is determined mainly by the diffusion process. 

Consider now the equation for lattice displacements (g^) with the force ([l9|). Note that in the hydrodynamic regime, 
Kv <C T~^, the dominant contribution to the force Gi comes from the local equilibrium partition function, i.e., the 
first term in (y), if we consider times greater than the electron-electron relaxation time, f 3> t. In this case, the force 
has the simple expression 



G,(r,i)=A, 



dxj 



where the constants 






are of the order of the electron density of states at the Fermi surface. 

From Eq. (Efl) with the help of the expression (|2l|) one can find the lattice deformation 
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(25) 



The first term in the brackets in Eq. (Ea) represents the particular solution of the inhomogeneous Eq. (M) while the 
second corresponds to the general solution of the homogeneous form of Eq. (|l|), and represents the effect of the surface. 
The integrand in (E5h contains poles associated with the diffuson and sound-wave excitations. Sound singularities are 
bypassed using infinitesimal phonon damping, uj —* uj + iO. 



V. EFFECT OF ACOUSTIC AND OPTICAL DISPLACEMENTS ON DESTRUCTION OF METALS 

Equation (E5^ describes the effect of nonequilibrium electron heating on lattice deformations of acoustic type. This 
deformation vanishes at the surface z = according to the boundary condition (O). For z ^ 0, the second term in 
brackets in ( |25| ) represents a deformation wave propagating from the surface into the bulk of the metal. It makes 
a nonzero contribution only at sufficiently small depths z < st ^ 10~^cm. Thus, we see that the deformation (Eq) 
peaks at z '--^ 10~^cm <C k~^. To obtain the order of the effect, we can drop the second term in parentheses. It is 
then convenient to integrate over lo, substituting the Fourier transform I{lu). We obtain 



duz 
dz 



Azzk{1~R) 

pP 



dt'I{t') 



dk 



-isk{t — t') 
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k^D^ + s2 



Akz 



(26) 



Consider times greater than the duration of a pulse t > to but less than the characteristic time of electron diffusion 
{t < (k^D)^^ ^ lO^^^s) and a sound-wave period {t <C (sk)^^ ^ lO^^^s) with characteristic wave vector of the order 
of the inverse skin depth k. In this range we can expand the exponentials in (Eq) in powers of t up to second order: 



duz 
dz 



AItoK{l - R) 2 



2p(3 



f f ^u{k)k^e'^ 
J 27r 
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Using the estimate A/p ~ gs'^/^'p and the Eq. 
result up to electron diffusion times t ~ (k^D)^^ : 



(|2J) for the maximum electron temperature, we extrapolate our 



dUz ,^ ^. Itp 

dz KJ3 



Tv^ep 



SI rr 



KTV^Ep 



(28) 



Setting s/v ^ 10 



10 cm , we arrive at the numerical estimate duz/dz ^ 10 g {Te / e f) 



Our result contains the natural factor T'^/e'j^, which means that laser heating is important as soon as the electron 
temperature is higher than the Fermi energy. Although the estimate was obtained for Te <SC £_f, it is still roughly 
correct up to T^ ^ Sf- The additional small factor s^ /{vkIY is due to the fact that the characteristic period of the 



sound wave {10~^^s) is much greater than the charachteristic times of electron diffusion {I0~^'^s) and laser heating 
(10~^^s). Therefore it would be of considerable interest to calculate the lattice deformation from high-frequency 
(but long-wavelength) excitations, i.e., optical phonons whose period is about IQ^^^s. This case differs from the 
calculations above in the equation of lattice motion (1) and electron force (2), due to the different form of deformation 
potential ( seecil ). Estimates show that the relative optical displacement (with respect to the lattice constant) are of 
order T^/ej,. 

VI. CONCLUSIONS 

Our result for acoustic deformation (Eq) agrees with the experiment reported by P. Rentzepis [17] where a defor- 
mation duz/dz ^ lO^'^ had been observed in the laser heating of noble metals. However, we see that the interaction 
of heated electrons with optical phonons can provide a more effective means of strong lattice deformation, but this 
case has yet to be studed experimentally. An ultrashort intense laser pulse can result in the destruction and ablation 
of metals, while only the electron component is heated, and the lattice stays at considerably low temperature. 

In conclusion, we emphasize two points. First, as follows from Eq. (9), the driving force for lattice expansion 
is proportional to TedTe/dz. Because of the high absorption coefhcient of metals in the UV (k ^ lO^cm^^), the 
temperature gradient reaches ^ 10^ K/cm. Note that the extremely high values of this parameter (which is typical 
of metals) leads to nonequilibrium expansion of the lattice. Second, subpicosecond elastic deformation of the lattice 
of order 10""^ — 10^^, corresponding to an internal pressure 10 — 100 GPa, can provide an effective mechanism for 
subsequent laser fracture of metals. 
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